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Raising Tc in charge density wave superconductor ZrTes by Ni 
intercalation 
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Abstract. - We report discovery of bulk superconductivity in Nio.osZrTea at Tc = 3.1 K, obtained 
through Ni intercalation. Superconductivity coexists with charge density wave (CDW) state with 
TcDW = 41 K. When compared to parent material ZrTes, filamentary superconducting transition 
is substantially increased whereas Tcdw was suppressed. The analysis of superconducting state 
indicates that Nio.osZrTea is an intermediately coupled superconductor. 
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Introduction. — Charge density vifave (CDW) and 
superconductivity are distinctive quantum orders that 
emerge due to Fermi surface (FS) instabilities |T]. Com- 
petition and coexistence between CDW and superconduc- 
tivity in low-dimensional materials is one of the most fun- 
damental problems in condensed matter physics. Among 
CDW superconductors, layered MX2 compounds and 
chain-typed MX3 compounds (M = transition metal, and 
X = S, Se or Te) are well-known [5]— [3]. It was recently 
pointed out that chalcogenide superconductors represent 
a weak coupling limit of the melting of the stripe (or smec- 
tic) order in cuprates [5]. Furthermore, the domelike Tc(x) 
[4] in Cua;TiSe2 and pairing mechanism was proposed to 
arise from the quantum criticality related to fluctuations 
in CDW order [6] — • Therefore discovery of new super- 
conductors in this materials class is of significant interest. 

ZrTes belongs to a family of trichalcogenides MX3 (M 
= Ti, Zr, Hf, U, Th, and X = S, Se, Te). The struc- 
ture consists of infinite X-X chains formed by stacking 
MX3 prisms along the crystallographic b axis. The poly- 
hedra are arranged in double sheets and stacked along 
monoclinic c axis by van der Waals forces (fig. 1(a)) 
[8]. ZrTes exhibits CDW transition at Tcdw = 63 K. 
The CDW nesting vector q cdw = (1/14, 0, 1/3) has no 
component in the chain-axis direction [9]. This is differ- 
ent from another well-known quasi-one-dimensional CDW 
compound NbScs with X-X chains which exhibits a re- 
sistivity anomaly along the chain axis [10]. The p(T) of 
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ZrTes is quasi two dimensional, metallic below 300 K with 
anomalies due to CDW transition [TT], and superconduc- 
tivity below Tc = 2 K. The superconductivity is not bulk 
but filamentary [T^]. Theoretical calculations [13] and 
photoemission study [T3] indicate that the CDW transi- 
tion is driven by the nesting of parallel planar sections of 
the FS which are related to the Te ^Px band along the Te- 
Te chains. Other parts of the FS remain unaffected and 
are responsible for superconductivity. Therefore, CDW 
and filamentary superconductivity coexist in ZrTes due 
to multiband nature of the Fermi surface. 

On the other hand, ZrTes (fig- 1(9')) exhibits uncommon 
pressure dependence of CDW and superconductivity [TS] . 
With increasing pressure Tcdw initially increases, then 
decreases up to 2 GPa and abruptly vanishes near 5 GPa. 
The Tc falls to below 1.2 K at 0.5 GPa and can not be ob- 
served up to 5 GPa. At higher pressures superconducting 
transition emerges again and increases continuously from 
2.5 K to 4.7 K at 11 GPa. The Tcdw{P) and Tc{P) above 
5 GPa can be explained by pressure induced FS modifica- 
tions. High pressures favor three-dimensionality, therefore 
the area of the planar portions of FS will decrease and the 
CDW transition will be suppressed. Consequently super- 
conducting Tc, a competing FS instability will be favored. 
However, the pressure dependence above 2 Gpa is still un- 
clear. The Tc(P) and Tcdw{P) indicate that the changes 
in FS are rather sensitive to crystal structure modifica- 
tion. Besides hydrostatic pressure, doping and interca- 
lation can be used to tune the crystal structure, lattice 
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properties, shape and characteristics of the FS. For exam- 
ple, intercalating external atoms into the inter layer weak 
van der Waals gap of MX2 can depress the CDW transi- 
tion and enhance the superconductivity [1]. Since there 
is also a van der Waals gap in ZrTea, intercalation is a 
powerful method to examine the relation between CDW 
and superconductivity in this system. Ni is often a suit- 
able intercalant as seen on the example of Nia;TaS2 and 
Ni:,TaSe2 [E], [II]. 

In this work, we report synthesis and physical prop- 
erties of Nio.osZrTes single crystals. Ni intercalation in- 
duces bulk superconductivity with enhanced Tc while the 
TcDW shifts to lower temperatures, indicating coexistence 
of CDW and superconductivity. 

Experimental. — Single crystal Nij-ZrTea was grown 
via the chemical vapor transport (CVT) method. The 
source and growth zone were set at 700 °C for 2 days 
and then kept at 720 °C and 645 °C respectively for 10 
days. Golden plate-like single crystals with a typical size 
of 2x2x0.2 mm'^ were obtained. The crystal structure 
and phase purity were examined by powder and single 
crystal X-ray diffraction pattern (XRD) with Cu Kq, ra- 
diation (A — 1.5418 A) using a Rigaku Miniflex X-ray 
machine. The structure parameters are extracted by fit- 
ting the XRD spectra using the Rietica software [18] . The 
composition of NiajZrTca single crystal was determined 
by examination of multiple points on the crystals using 
energy dispersive X-ray spectroscopy (EDS) in an JEOL 
JSM-6500 scanning electron microscope. The measured 
compositions are Nip. o52(3)Zrioo(2)Te3. 10(2), abbreviated 
as Nio.osZrTca. Physical property measurements were per- 
formed in a Quantum Design PPMS-9 and MPMS XL 5 
instruments. 

Results and discussion. — Powder X-ray diffraction 
(XRD) result (fig. 1(b)) of ground samples indicates that 
all peaks of Nio.osZrTes can be indexed using the struc- 
ture of ZrTea, i.e. the intercalation does not change the 
structure of the mother compound. The fitted lattice pa- 
rameters are a = 0.5895(1) nm, b = 0.3923(1) nm, c — 
1.0160(2) nm, and /3 = 97.76(1)°. Our refinement of pure 
ZrTca crystals gave lattice parameters a ~ 0.586(3) nm, 
b = 0.392(7) nm, c = 1.009(5) nm and /? = 97.75(1)°. 
Hence, the a and c - axis lattice parameters of Nio.osZrTcs 
are larger than those of ZrTca . The substitution on Zr site 
is unlikely due to the smaller size of Ni ion (r7Vi2+ — 69 
pm) than Zr iow[v zr'^+ = 72 pm). This indicates that Ni is 
most likely intercalated in the Van-der- Waals bonded crys- 
tallographic layers in ZrTca. The XRD pattern of single 
crystals (ig. 1(b) inset) reveals that the crystal surface is 
normal to the c axis with the plate-shaped surface parallel 
to the a6-plane (fig. 1(c)). 

Figure 2(a) shows the temperature dependence of resis- 
tivity at zero field from 1.9 K to 300 K. Both Pa(T), and 
Pb(T) undergo a relatively sharp superconducting transi- 
tion at Tconset = 3.74(3) K and 3.59(2) K, respectively 
(upper inset). It should be noted that the Pb{^) is metal- 
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Fig. 1: (a) The structure of ZrTea with green (orange) symbols 
for Zr (Te) atoms, (b) Powder XRD pattern of Nio.osZrTea, 
Inset: photo of typical single crystals of Nio.osZrTea. (c) XRD 
result of Nio.osZrTes single crystal. 




Fig. 2: (a) Temperature dependence of the resistivity paiT), 
PbiT) of Nio.osZrTea. Inset: resistivity near Tc. (b) Temper- 
ature dependence of ac magnetic susceptibility of Nio. 05 ZrTea. 
Inset: temperature dependence of dc magnetic susceptibility 
with ZFC and FC. 
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lie in the normal state with no anomalies due to CDW and 
with a residual resistivity ratio (RRR) of 6. The /9a(T) ex- 
hibits a transition located at 41 K (lower inset fig. 2(a)). 
This can be ascribed to CDW since it is similar to resis- 
tivity anomalies in ZrTea [TT]. However, the transition 
has been depressed from around 65 K to 41 K. Therefore, 
there is a coexistence of superconductivity and CDW in 
Ni2:ZrTe3 similar to ZrTca but with nearly double bulk 
Tc [12]. The filamentary superconductivity and the CDW 
originate from different sections of FS in ZrTea. From 
the difference in T^ Tcdw in Nio.osZrTea and ZrTea we 
can conclude that Ni intercalation dramatically changes 
the FS related to the superconductivity but has a minor 
effects on FS parts associated with the CDW transition. 
Intercalation has a distinctively different effect when com- 
pared to pressure. The reason could be that intercalation 
only modifies one part of FS related to superconductiv- 
ity whereas pressure increases three-dimensionality of all 
parts. 

Figure 2(b) shows the temperature dependence of the ac 
susceptibility of Nio.osZrTes single crystal with dimension 
of 2x1.5x0.24 mm'^ for H||ab and H||c. The single sharp 
peak of 47rx" accompanied by a very steep transition in 
47rx', indicates that the sample is rather homogeneous. 
The onset temperature for both field direction is 3.1 K. 
This is slightly lower than that obtained from resistivity 
measurement. The transition width ATc is 0.3 K. For 
H||ab, the value of -47rx' at 1.8 K is 110% without any 
demagnetization factor correction, indicating that the su- 
perconductivity is bulk. This is different from the mother 
compound ZrTca where superconductivity is undetectable 
in magnetization due to its filamentary nature below T^ = 
2 K. Inset in fig. 2(b) shows the dc magnetic susceptibil- 
ity for field along a, c axes with zero-field cooling (ZFC) 
and field cooling (FC). For ZFC curves, the Tc^onset and 
ATc are consistent with the results of ac susceptibility and 
the values of -47rx at 1.8 K are nearly identical. On the 
other hand, the volume fraction estimated from the FC 
curve is about 2% at 1.8 K for II||c, similar to other inter- 
calated compounds such as (Pyridine) i/2TaS2 [19], YbCg 
[20] . Nia;TaS2 [21]. The small magnetization values for 
FC is likely due to the complicated magnetic flux pinning 
effects in the intercalated compounds [19] . 

The results of the dc magnetization versus field M(II) 
at various temperatures for both directions are shown in 
fig. 3(a). The shape of the M(H) curves points that 
Nio.osZrTca is a typical type-II superconductor. Figures 
3(b) and (c) show the initial M(H) curves at the low field 
region. All curves clearly fall on the same line and devi- 
ate from linearity for different temperatures. Linear fits 
describe the Meissner shielding effects ("Meissner line") 
(figs. 3(b) and (c)) The obtained slope of the linear fit up 
to 30 Oe at the lowest temperature of our measurement T 
= 1.8 K is -0.995(5)« -1. This corresponds to -47rM = H 
for H||ab where the demagnetization factor is negligible. 

The value of H*^ at which the field starts to penetrate 
into the sample can be determined by examining the point 
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Fig. 3: (a) Magnetization hysteresis loops of Nio.osZrTes at T 
= 1.8 K for H||ab and H||c. (b) and (c) Low field parts of M(H) 
at various temperature for H|[ab and H||c with demagnetization 
correct, respectively. The two solid red lines are the "Meissner 
line" as discussed in the text, (d) Temperature dependence 
of Hci for H||ab and H||c. The dashed lines are the fitted 
lines using H^i ~ J/ci(0)(l-(T/Tc-)^). Inset: the temperature 
dependence of anisotropy of Hd, 7^^ = Hci,ab{T) / Hci,c{T) . 



of deviation from the Meissner line on the initial slope of 
the magnetization curve. The H*^ is not the same as the 
real lower critical field, Hd, due to the geometric effect. 
The Hci can be deduced from the first penetration field 
H*i, assuming that the magnetization M = -Hd when 
the first vortex enters into the sample. Thus H has been 
rescaled to H = Ha-NM and Hd = H*^/{1-'N) where N 
is the demagnetization factor and Ha is the external field 
[22] . We estimate demagnetization factors 0.075 and 0.781 
for H||ab and H||c by using Hd = H*^/teiiih{y/0 Mb /a), 
where a and b are width and thickness of a plate-like su- 
perconductor [23]. This is consistent with previous, i.e. 
the demagnetization factor for H||ab is negligible. How- 
ever, demagnetization factor must be considered for H||c. 
The M(H) curve for H||c with demagnetization correction 
is shown in fig. 3(c) and the slope of the fitted line in 
linear region is -0.963(1) « -1. From the above and con- 
strained with the resolution limit of our magnetization 
measurement Am= 3x10^^ emu we extract the tempera- 
ture dependence oi Hd for both field directions (fig. 3(d)). 
The Hd{T) can be well explained by Hci(T) = i?ci(0)[l 
- {T/Tc)% The obtained i?ci,ah(0) and i/ci,c(0) are 50(1) 
and 321(8) Oe, respectively. On the other hand, the 
anisotropy of Hd, 7h,i = HdM//c/Hd,H//ab, increases 
with decreasing temperature and is estimated about 6.4 at 
T = 1.8 K. This is much larger than in MX2 compounds 

m- 

The main panel of fig. 4(a) shows the temperature de- 
pendence of the specific heat of Nio.osZrTea below 5 K at 
H = 0, 10, and 50 kOe. A jump at 3 K can be clearly seen 
at H = 0, indicating bulk superconductivity of Nio.osZrTea 
sample. The Tf, onset of 3.1 K deternnned from the specific 
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Fig. 4: (a) Low temperature specific iieat beiiavior of 
Nio.osZrTea plotted as Cp/T vs T at H = 0, 1 and 50 kOe. 
Tlie solid line is a fit described in the text. Inset: temperature 
dependence of Cp(T) from 1.95 K to 300 K at H = kOe. (b) 
Temperature dependence of the electronic specific heat plotted 
as Ces/T vs T at H = kOe. The solid line shows fitted result 
of Ces/T assuming an isotropic s-wave BCS gap. 



heat jump is consistent with that obtained from the sus- 
ceptibihty and transport measurements. At H = 1 kOe, 
the superconducting anomaly is shifted to lower temper- 
atures. The normal state is recovered at H = 50 kOe. 
In order to obtain the electronic specific heat coefficient 7 
and Debye temperatures Qd , the low temperature specific 
heat at H = 50 kOe is fitted using Cp/T = 7 + I3T^. The 
obtained parameters are 7 = 2.7(1) mJ/mol-K^ and Qd 
= 192.4(1) K using Qd = {12ti'^ NR/hfif/^ where N = 4 
is the number of atoms per formula unit and R is the gas 
constant. The high-temperature specific heat approaches 
the value of 3NR, in accordance with the Dulong-Petit 
law (fig. 4(a) inset). According to the McMillan formula 
for electron-phonon mediated superconductivity [53], the 
electron-phonon coupling constant A can be determined 

by 



Tr 
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1.45 



exp[ 



1.04(1 + A) 
A-Ai*(l + 0.62A)^ 



(1) 



where /i* ~ 0.13 is the common value for Coulomb pseu- 
dopotential. By using T^ = 3.1 K and Qd = 192.4 K, 
we obtain A = 0.63, indicating an intermediately coupled 
superconductor. The electronic specific heat part Ces ob- 
tained by subtracting phonon part from the total specific 
heat is shown in fig. 4(b). From fits to Ces/T-T below 
Tc using the BCS formula for the electronic contribution 
to the specific heat [25] , the ratio of the gap and the crit- 
ical temperature is about, 2A/A;bTc = 3.00(2). This is 
smaller than the typical BCS value (3.53) in the weak- 
coupling limit [3S]. The solid red line in fig. 4(b) shows 
the data simulation for an isotropic s-wave BCS gap. The 
specific heat jump at Tc, ACes/lTc — 1.20(3), is somewhat 
smaller than the weak coupling value 1.43 [24]. The cal- 
culated coupling strength from the heat capacity anomaly 
could be influenced by the rather large fitting inaccuracy 
near Tc and somewhat broad superconducting transition 
width in C/T{T). 

Temperature dependent resistivity of PbC^) of 
Nio.osZrTca below 5 K in magnetic fields for II||a 
and H|jc is shown in figs. 5(a) and (b). With increasing 




Fig. 5: Temperature dependence of the resistivity pb{T) of 
Nio.osZrTea for (a) H||a and (b) H||c at the various magnetic 
fields, (c) Temperature dependence of the upper critical field 
Hc2 for H||a and H||c. Solid lines show the fitted result using 
WHH formula described in the text. 



magnetic fields, the resistivity transition width becomes 
broader and the onset of superconductivity gradually 
shifts to lower temperatures. At H = 5 kOe, for H|jc, the 
superconducting transition can not be observed above 
1.9K, whereas for H||a, the superconductivity above 1.9 
K is nearly suppressed at H = 10 kOe. 

Figure 5(c) shows the upper critical field Hc2{T) of 
Nio.osZrTcs corresponding to the temperatures where the 
resistivity drops to 90% of the normal state resistivity 
Pn,b{T, H){Tc^onset) at 5 K. Since Pauli limiting field Hp{Q) 
= 1.84rc > 55 kOe, the orbital effect should be the domi- 
nant pair-breaking mechanism. According to the conven- 
tional one-band Werthamcr-Helfand-Hohenberg (WHH) 
theory, which describes the orbital limited upper critical 
field of dirty type-II superconductors pS] , the Hc2 can be 
described by 
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T /Tc, ■0 is a digamma function and 

^Hc2 



(2) 



(3) 



^^Tc{-dHc2/dT)T=T^' 

The Hc2{T) fits are shown by solid and dotted lines in 
fig. 5(c). The obtained -dHc2/dT\T=n = 5.049(285) 
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Table 1 


Superconducting parameters of Nio.osZrTes. 






Nio.osZrTea 


HcM (Oc) 


i/e2(0) (kOe) iJ,(0) (Oc) m (nm) 


A(0) (nm) 


k(0) 


H||ab 

H||c 


50(1) 
321(8) 


12.281(860) 449(13) 26(1) 
4.862(340) 10(1) 


78(11) 
1242(100) 


19.3(8) 
3.0(3) 



kOe, Tc = 3.51(4) K for H||a and -dHc2/dT\T=T, = 
2.016(122) kOe, T^ = 3.48(4) K for H||c are consistent 
with the Tconset [pb = 90%p„,6(T, H = 0) = 3.59(2) 
K). The estimated upper critical fields are iJc2(0) = 
12.281(860) kOe and 4.862(340) kOe for H||a and H||c. 
From the i?c2(0) zero-temperature coherence length ^(0) 
can be estimated with Ginzburg-Landau formula i?c2,c(0) 
= $o/[2<^,(0)], and iJc2,a6(0) - $o/[2^eab(0)ec(0)] 
where $o — 2.07x10^^^ Wb. Based on the values of 
i?ci(0) and ifc2(0), the Ginzberg-Landau (GL) parame- 
ter K is obtained from iJc2(0)/i?ci(0) = 2k'^ /{biK + QM). 
And thermodynamic critical field iJc(O) can be obtained 
from iJc(O) = i?ci,Qf,(0)/[%/2Kab(0)]. The GL penetration 
length Aab(O) and Ac(0) can be evaluated using Kc(0) = 
Aab(0)/eab(0), and Kab{Q) = [Aa6(0) Ac(0)/U(0)6(0)] '/' 
[27] . All of obtained parameters are listed in Table 1. 

Conclusion. — In summary, we have discovered the 
bulk superconductivity in Ni intercalated ZrTea and pre- 
sented detailed characterization of superconducting state. 
The bulk superconducting Tc was nearly doubled when 
compared to filamentary Tc in parent material ZrTcs. The 
normal state electrical resistivity indicates that the CDW 
transition is somewhat suppressed with Ni intercalation, 
coexisting with bulk superconductivity. The Ni intercala- 
tion substantially changes the FS related to the supercon- 
ductivity but has a minor effects on FS parts associated 
with CDW. This is rather different from the pressure ef- 
fects. Our results suggest intermediately coupled super- 
conductivity and possible rich vortex physics. 
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